Spin-3/2 random quantum antiferromagnetic chains 
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We use a modified perturbative renormafization group approach to study the random quantum 
antiferromagnetic spin-3/2 chain. We find that in the case of rectangular distributions there is a 
quantum Griffiths phase and we obtain the dynamical critical exponent Z as a function of disorder. 
Only in the case of extreme disorder, characterized by a power law distribution of exchange couplings, 
we find evidence that a random singlet phase could be reached. We discuss the differences between 
our results and those obtained by other approaches. 



The study of random antiferromagnetic chains is an 
important and actual area in magnetism. Since by now, 
many of the physical properties of pure chains are under- 
stood, it is natural to include disorder in these systems 
and look for the modifications it introduces. In the case of 
spin- 1/2 random exchange Heisenberg antiferromagnetic 
chains (REHAC) a perturbative approach introduced by 
Ma, Dasgupta and Hu (MDH) was very successful 1] to 
investigate these systems. This approach turns out to 
be asymptotically exact and this allowed Fisher ^ to 
fully characterize the properties of the new disordered 
phase, for which, the name random singlet phase was 
coined. Unfortunately when generalized to higher spins 
this method, in its simplest version at least, revealed 
to be ineffective. The reason is that in the elimination 
procedure of the strongest bond Q, the new interaction 
between the spins, coupled by exchanges Ji and J2 to 
the strongest coupled pair is given by, J' = {2/3)S{S -\- 
l)JiJ2/^ I- For 5 > 1 the factor (2/3)5(5' 1) > 1 
and the problem becomes essentially non-perturbative for 
arbitrary distributions of exchange interactions. Several 
approaches have been introduced to circumvent this diffi- 
culty 4, 5,^I.,^,l|il|ll[li[ll[l4, 15,Ifi.] which how- 
ever not always lead to the same unambiguous results. In 
this Communication we apply a previous method used to 
treat the random spin-1 chain (13.] for the spin-3/2 RE- 
HAC. This particular chain has been the subject of recent 
studies 0, 0, 01 and it would be nice to confirm the 
results obtained in these works using another approach. 



As mentioned previously the method of Ma, Dasgupta 
and Hu consists in finding the strongest interaction (Q) 
between pairs of spins in the chain (see Fig. la) and treat- 
ing the couplings of this pair with its neighbors ( Ji and 
J2) as a perturbation. For a chain of spins S = 3/2, 
after elimination of the strongest coupled pair, the new 
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FIG. 1: The two elimination procedures as described in the 
text (Jl > J2). 



coupling between its neighbors is given by 



, _ 5 J1J2 



(1) 



Consider the case Ji > J2. If Ji > (2/5)ri than the 
new effective interaction J' is necessarily larger than one 
of those eliminated, in this case, than J2. 

Our generalization of the MDH method consists in 
either of the following procedures shown in Fig. ^ If 
the largest neighboring interaction to Q,, Ji < (2/5)57, 
then we eliminate the strongest coupled pair obtaining 
an effective interaction between the neighbors to this 
pair which is given by Eq. ^ (see Fig. la). This new 
effective interaction is always smaller than those elim- 
inated. Now suppose Jl > J2 and Ji > (2/5)17. In 
this case, we consider the trio of spins 5 = 3/2 cou- 
pled by the two strongest interactions of the trio, Ji and 
f2 and solve it exactly (see Fig. lb). The ground state 
of this trio of spins 5 = 3/2 is a degenerate quadru- 
plet. It will be substituted by an effective spin— 3/2 in- 
teracting with its neighbors through new renormalized 
interactions obtained by degenerate perturbation theory 
acting on the ground state of the trio. This procedure 
which implies diagonalizing the 64X64 matrix of the trio 
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FIG. 2: The evolution of the cut-off, for an initial rectangular 
distribution with G — 0.5, along the renormalization process 
of the spin-3/2 REHAC. We show the results for the origi- 
nal MDH and the present (modified) renormalization group 
procedures. 



is carried out analytically. This is important for obtain- 
ing results on large chains and to deal with the large 
numbers of initial configurations that we use. These 
procedures guarantee that we always comply with the 
criterion of validity of perturbation theory as shown in 
Fig. 121 We have considered initial rectangular distribu- 
tions, Pa{J) = - G))e{n - J)e(J - G) of inter- 
actions and even for the weak disorder case, with a gap 
G as large as G = 0.5 in the distribution (fi = 1), the 
method works very well and never an interaction larger 
than those eliminated is generated (see Fig.EJ. 

The phase diagram for rectangular original distribu- 
tions can now be obtained. For strong disorder, which 
corresponds to G = 0, we find a Griffiths phase with a 
dynamic exponent Z ^ 12.7 as shown in Fig. 13 This 
phase is characterized by first gap distributions that sat- 
urate at low energies in the form, P(— log A) ^ t^l^ for 
This is obtained starting from a given 
configuration of random interactions for a chain of size 
h and eliminating the spins, as described above, until a 
single pair remains. The interaction between these re- 
maining spins yields the first gap A for excitation. The 
dynamic exponent Z relates the scales of length and en- 
ergy through A oc . 

In Figure 121 we show the first gap distributions for dif- 
ferent degrees of disorder as characterized by different 
gaps G in the initial distribution of interactions. For all 
cases, including that of strong disorder (G = 0), we find 
that the first gap distributions saturate at low energies 
with the dynamic exponent Z independent of L for L 
sufficiently large. We have to consider large chains in 
order to observe this effect. We find Z^ ~ 0.43 and 
Zoo ~ 1.12 for G = 0.2 and G = 0.12, respectively. From 
these values of the dynamic exponent we can deduce the 
existence of a Griffiths phase extending up to Gc ~ 0.11 
where the dynamic exponent reaches the value Z = \. 
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FIG. 3: Probability distributions of first gap for initial rect- 
angular distributions of couplings with gaps G and different 
systems sizes L. For clarity not all values of L are shown. The 
solid lines represent best fits to the form logj^Q[P(— logj^Q A)] = 

= 12.68, 



Al- logio A- a) G = 0, Z3, 



Z]f i"SlO a.) ^ — u, ^3000 = 10.51, Z700O 

and Zgooo = 12.70. b) G = 0.12, Z3000 = 0.87, Z7000 = l-H 
and Zsooo = 1.12. c) G = 0.2, Z3000 = 0.35, ^'tooo = 0.43 
and ^8000 = 0.43. 
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Since, for example, the susceptibility x ^ T^^^ , a singu- 
lar low temperature behavior implies Z > 1. At Gc there 
is in fact a significant change in the nature of the ther- 
modynamic behavior of the system jl6| . The phase for 
G > Gc is one with quasi- long range order, i.e., with spin 
correlations decaying algebraically with distance, similar 
to the zero temperature phase of the pure chain [l^. 

In order to check the existence of a random singlet 
phase in the spin-3/2 chain we consider another class of 
distributions of exchan ge c ouplings associated with ex- 
treme disorder 12(1 . l2ll . These distributions are of 
the form, P(J) oc J~^^^/^. For 5=1 this reduces to 
the gapless case of rectangular distributions considered 
previously and for S > 1 , we have the extreme disordered 
cases. We now report our results for the random spin-3/2 
chain obtained with the modified renormalization group 
procedure for the case of an extreme disordered dis- 
tribution with 6 = 20. In a random singlet phase the 
fixed point distribution of interactions which is attained 
when the cut-off ft is sufficiently reduced, takes the form 



(2) 



The exponent a is a function of the cut-off Q and varies 
as, a = — l/lnfi. Also for a random singlet phase the 
fraction of remaining active spins o as a function of the 
energy scale set by the cut-off [3 is given by. 
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FIG. 4: a) The exponent a describing the asymptotic low 
energy behavior of the renormalized exchange distribution as 
a function of the scale dependent cut-off. b) Fraction of active 
spins as a function of the cut-off. For comparison we show the 
results for the spin-1 and spin-3/2 REHACs. In both cases 
the starting distribution is extremely disordered with 5 = 20 
(see text). 



with spin- 1/2 degrees of freedom. These results are quite 
distinct from those obtained above where a random sin- 
glet phase is hardly evident even for extremely disordered 
original distributions. 

The decomposition of a chain of spins-S in smaller 
spins, relies on projecting out the highest energy level of 
a pair of spins-S*. However the remaining excited states 
are kept in this procedure to maintain the correct num- 
ber of states. For example, in the case of a pair of spins-1 
with a total of nine states, the singlet ground state and 
the first excited triplet are kept to yield the four states of 
the relevant antiferromagnetically coupled spin-1/2 pair 



^ L llnf^li/V-- 



(3) 



The exponent ■0 establishes the connection between the 
characteristic length L and the energy scale fi. This is an 
extension of the usual definition of a dynamic exponent 
{Qr^ (X T (X L^) for the case of logarithmic scaling [2^ . 
In Fig. 2Ji we show the exponents a obtained from the 
asymptotic form of the exchange distributions after the 
cut-off n has been sufficiently reduced. For comparison 
we show the results for a spin-1 REHAC with the same 
original extreme disorder distribution for which a random 
singlet phase is clearly established. We have considered 
here chains of size as large as L = 4.5 x 10^. In Fig.^ we 
show the density p = 1/L of active spins as a function 
of the cut-off ri. From this expression we extract the 
exponent ip (see EqEJ which takes the value V' = 1/(2-4) 
close to the value -0 = 1/2 expected for a random singlet 
phase 0. As shown in this figure, for comparison, the 
spin-1 chain has clearly converged to this phase within 
the same scale of the cut-offs. Our results suggest that 
in this case of extreme disorder, the spin-3/2 REHAC 
eventually reaches a random singlet phase, although the 
convergence is very slow. 

Recently another approach to the spin-3/2 chains has 
predicted the existence of a random singlet phase in these 
chains, even for weak disorder [l3 . fl5i |. This is associated 



The MDH elimination procedure can be generalized for 
finite temperatures and arbitrary spins S. It is given by 



J' = -SiS- 



(4) 



where 



Ws{y) = 



{2S+ 1)2 -E-lf (2*+ l)e-5H''+i)2' +1)] 



45(5 + l)E:=f(2* + l)e-^ 



(5) 

Notice that for sufficiently high temperatures, the factor 
p{S + l)Wsif3n) < 1 and the MDH elimination pro- 
cedure works in this case. A random singlet phase is 
reached in the sense that the asymptotic distribution of 
exchange attains the form given by Eq.|21at these temper- 
atures. However, as T is reduced the problem becomes 
essentially non-perturbative, for spins S > 1, as the equa- 
tion above generates coupling larger than those elimi- 
nated. In particular at T = the excited states which 
reduce the factor Ws from its value Ws{T = 0) = 1 
are now frozen. In fact none of the excited states play a 
role in the problem at zero temperature. Notice that in 
our generalized renormalization scheme, degenerate per- 
turbation theory is applied to the ground state of the 
spin trio. We believe this is the main reason for the dis- 
crepancy between our results and those obtained by the 
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FIG. 5: Temperature T* below which the MDH perturbation 
theory breaks down for different values of the spin S. The en- 
ergy fesT* is in units of the cut-off SI of the original exchange 
distribution. 



glet phase, associated with spin-1/2 degrees of freedom, 
even for weak disorder. We attribute the difference be- 
tween these results and those we have obtained to the 
fact that the former approach takes into account excited 
states which mimic the effects of temperature and favor 
the appearance of a random singlet phase. Our results 
however are consistent with those of Carlon et al. 0| 
that find a random singlet phase in spin-3/2 chains for 
the case of extremely disordered distributions. 

We would like to thank Conselho Nacional de 
Desenvolvimento Cientifico e Tecnologico-CNPq-Brasil 
(PRONEX98/MCT-CNPq-0364.00/00), Fundacao de 
Amparo a Pesquisa do Estado do Rio de Janeiro- 
FAPERJ for partial financial support. 



authors of Refs |14L Il5| . The consideration of excited 
states in the problem favors the appearance of an infi- 
nite disorder random singlet phase, as occurs at finite 
temperatures. 

In the limit 5 — > cxd and T — > 0, replacing the sum 
by an integral and with a proper renormalization of the 
Hamiltonian, Eqs. 0] and |S1 yield J' = J1J2/ '^ksT , in 
agreement with the result of Ref Q for classical spins. 
In Fig. |S1 we show the temperature T* below which the 
simple perturbative approach breaks down, for a given 
value of the spin of the random chain. 

Notice that for rectangular distributions, the Griffiths 
phase of the spin-1 REHAC extends up to Gc = 0.45 and 
for spin-3/2 up to Gc — 0.11. For random classical spin 
chains, the susceptibility x ^ ^(0)|lnT| where P{Q) is 
the finite weight at the origin of the original distribution 
P{J) 0. This weak logarithmic singularity is similar to 
that expected for quantum chains at the border of the 
Griffiths phase {Z k, 1) as if in this case, Gc = 0. 

We have studied a spin-3/2 REHAC using an exten- 
sion of the renormalization group procedure introduced 
by Ma, Dasgupta and Hu 1]. This method which con- 
siders larger clusters of spins eliminates the difficulties 
associated with the perturbative nature of the MDH pro- 
cedure for the cases of spins S >\. The new procedure 
works very well for the spin-3/2 chain and never bonds 
larger than those eliminated are generated for the dis- 
tributions used here. For rectangular distributions we 
found the spin-3/2 REHAC presents a Griffiths phase 
up to a critical value of disorder. We have also con- 
sidered the case of extreme disorder, where the starting 
exchange distributions are singular for small values of the 
coupling. For values of the disorder parameter as large 
as (5 = 20 our results only suggest that a random singlet 
phase will be asymptotically reached as the cut-off of the 
distribution f2 — > 0. We have compared our results with 
those of another approach which predicts a random sin- 
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